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ABSTRACT 



The propulsion losses, which are caused by added drag 
due to steering of the ship, can be minimized by using an 
Adaptive Automatic Optimal Controller. It is shown in this 
thesis that an Adaptive Automatic Optimal Controller is 
capable of providing fuel savings in excess of 0.5 % over a 
well tuned PID controller when operating at the design speed 
at random headings in sea states. A new approach was used 
in finding fuel savings without using the engine specifica- 
tions . 

It is shown that the second- order forces and moments 
create drift motions along the surge, sway, and yaw axes. 
As a consequence of this, the second-order forces and 
moments cause more fuel consumption than the first-order 
forces and moments, which create only oscillatory ship 
motions along these axes. So the sea state in the determin- 
istic model is represented by the first-order and second- 
order forces and moments. 
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I. INTRODUCTION 



An overall rise in fuel prices has led to an increasing 
interest in the design of autopilots for ships. The purpose 
of the automatic steering control is to minimize propulsion 
losses which are caused by added drag due to steering of the 
ship . 

Several attempts have been made to define a measure of 
steering efficiency based on propulsion losses and proposals 
have been made for the design of an autopilot which mini- 
mizes the propulsion losses. As a performance criteria, 
added resistance due to steering suggested by Norrbin, N. H. 
was used in most of the studies because it is convenient for 
ship board use. 

Many researchers indicate that a carefully designed 
controller could save from one to two percent of fuel by 
minimizing the propulsion losses which are caused by added 
drag due to steering of the ship . No analytical evaluation 
of propulsion losses due to steering in a sea state has been 
made . 

The goal of this thesis is to study fuel savings using 
various kinds of controllers, and especially to compare an 
Adaptive Automatic Optimal Controller with a PID controller 
to see the difference between them. The optimal gain parame- 
ters of the Adaptive Automatic Optimal Controller are 
provided by LTJG. Cetin Diken [Ref. 1]. 

To study the optimization problem, models of both the 
ship and the operating environment are required. Chapter 2 
addresses what type of computer model can be used to repre- 
sent the ship. 

The ship's nonlinear equations of motion were needed to 
simulate the ship in the computer program. Chapter 3 
addresses the Mariner Class ship nonlinear equations of 
motion . 
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Chapter 4 addresses the problems of testing the ship 
simulation model in calm water, of expressing the forces 
and moments due to sea state, and the effects of sea state 
on ship's behaviour. 

Chapter 5 addresses the derivation of the cost function 
which represents the added drag due to steering. 

What is an adequate way to represent the fuel consump- 
tion ratio of added drag due to steering? Chapter 6 
addresses the fuel consumption ratio. 

Ship dynamics change with operating conditions such as 
ship speed, encounter angle, encounter frequency, and sea 
state. Chapter 7 presents the fuel consumption ratio of 
added drag due to steering for PID and an Automatic Optimal 
Controller, and the difference between them. 

Conclusions were drawn from simulation results, and are 
presented in Chapter 8. This chapter also recommends topics 
for future studies. 
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II. DESCRIPTION OF THE COMPUTER MODEL 



Before a specific controller can be designed, a real- 
istic model for the ship dynamics must be found to enable 
proper simulation of the ship movements in response to 
control signals. 

The model which best represents ship- steering dynamics 
is a Taylor's series expansion of the force and moment rela- 
tionship around a selected steady-state operating point. The 
resulting equations are commonly known as the equations of 
motion .[ Ref s . 2, 3, 4, 5] 

A computer program was developed in order to provide a 
computer simulation for the ship, using experimentally meas- 
ured hydrodynamic coefficients for the Mariner Class Ship 
[Refs. 2, 3, 6, 7]. Figure 2.1 shows the block diagram. The 
computer program is shown in Appendix C [Refs. 8, 9, 10, 
11 ] . 

The function minimization subroutine used was BOXPLX 
which was programmed by R. Hilleary. The task of BOXPLX is 
to find the minimum of any function. It may be subjected to 
explicit constraints of the variable or implicit constraints 
on functions of the variables. [Refs. 11] 

The sea disturbance is found as first and second-order 
forces and moments by a sea state program which has been 
written by J. Cass [Ref. 9]. 

D type or PID controller was used in the computer simu- 
lation program. PID controlller is shown in Appendix A. The 
D type controller was described by: 



<5 = 



K(l+T z s) 



d + T p s) 



T I S 



*e 
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where K,T z ,Tp, and Tj are optimal parameters which were 
found by the minimization subroutine, BOXPLX. The defini- 
tions of symbols used in equations and figures are given in 
Appendix D. 




Figure 2 . 1 



Block Diagram of Ship and Control System. 



Ill . NONLINEAR EQUATIONS OF MOTION 



Linear theory is useful for analyzing the influence of 
ship features on controls-fixed stability as well as on the 
turning ability of stable ships in the linear range [Refs. 
2, 3, 4, 5, 12, 13 ] . However, it fails to predict accurately 

the characteristics of the tight maneuvers that most ships 
are capable of performing and it cannot predict the maneu- 
vers of unstable ships. 

Nonlinear equations of motion are Suitable for 
predicting tight maneuvers and also suitable for computer 
programming. The nonlinear equations of motion based on a. 

'Taylor series expansion of forces and moments including 
terms up to the third order have been developed by Abkowitz 
[Refs. 5, 12], and Strom-Tejsen [Ref. 3]. Accuracy is not 

improved by including terms higher than third order. The 
development of these nonlinear equations is based on a 
restatement of linear equations to include rudder angle. 
Equations X, Y, and N are functions of u, v, r, u, v, r, and 

5 . 

Combining the nonlinear Taylor series expansion of 
forces and moments terms up to third order with the dynamic 
response terms of the X, Y, and N equations, the nonlinear 
equations of motion are shown in Appendix B. The hydrody- 
namic coefficients of the equations are determined from 
experimental data obtained from captive model tests. They are 
given in Appendix B [Ref. 7]. The terms not included in 
tables are negligible. 
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IV. SHIP'S BEHAVIOUR IN CALM WATER AND SEA STATES 



A . CALM WATER 

A simulation program in Appendix C was run for turning 
and zig-zag maneuvers to observe u, v, r, and \p using the 
ship's nonlinear equations of motion and Mariner Class Ship 
coefficients . 

It was observed that the rudder angle changes the ship's 
course. The ship's speed decreased while turning. The abso- 
lute value of v and r increased and after 'a while they 
reached steady- state values, since the ordered rudder angle 
is constant. The larger the rudder angle, the greater the 
decrease in speed, and the faster the ship turns. 

Time responses of r, v, and u are shown in Figure 4.1, 
4.2, 4.3, and X vs . Y is shown in Figure 4.4 for the turning 
maneuver with 5 = 25 degree. S, , and u are shown in Figure 
4.5 and Figure 4.6 for the zig-zag maneuver. As is seen from 
the figures, a suitable and sufficiently accurate ship 
computer model was defined. 

B . SEA STATE 

To observe the ship's behavior in a sea state, distur- 
bance forces and moments are needed. They depend on sea 
state, ship speed, encounter angle, and encounter frequency. 
The added mass and added inertia are functions of encounter 
frequency and sea state. 

A regular sea model was used as the sea representation: 
the wave crests assumed to be straight, infinitely long, 
parallel, and equally spaced with constant wave height 
[Refs. 2, 9, 14, 15]. 

The forces exerted by the regular sea can be represented 
as the sum of two components, called first and second-order 
forces [Refs. 2, 6, 9, 14, 16, 17, 18, 19, 20]. 
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The first-order exciting forces are a zero mean process. 
The second-order forces are a nonzero mean process, i.e. , 
they have essentially constant values. In regard to maneu- 
vering, the first-order forces are not the primary mover of 
the ship. The reason for this is that the high frequencies 
of the first-order forces are higher than the ship can 
readily respond to. The second-order forces cause the large 
excursions that must be manually controlled, while the 
first-order forces cause only a ripple on the ship surface 
trajectory [Ref. 6]. 

The sea also induces moments on the hull of a surface 
ship. These can be represented by first and second-order 
moments which possess the same basic characteristics as the 
first and second-order forces, and contribute to angular 
motion in a similar way. 

The first-order forces exerted by the regular sea have 
the form: 



F = w a R i Cos(<o e t+<£ i ) 



( eqn 4.1) 



where 



2 . 

— ^wave /*?) [ ^/^wave) ” u Cos^J) ] 

The second-order forces FXX, FYY, and moment MZZ have 
constant components and periodic components at twice the 
encounter frequency [Refs. 6, 9, 20]. The second-order 

forces FXX, FYY, and moment MZZ were assumed to approximate 
the constant drifting forces and moments calculated by using 
Sea State Program of James L. Cass [Ref. 9] for the Mariner 
Class ship in regular waves. This approximation was- used 
because these were the data available at the time. The 
sinusoidal part of the the second-order forces and moments 
(frequency of 2 w e ) was neglected at this point of the work 
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[Ref. 6]. The high frequencies of the second-order forces 
and moments are also’ higher than the ship can respond to. 
They do not effect the ship surface trajectory. 

The exciting forces Rj_ , and second-order forces for 
different encounter frequencies and encounter angle were 
obtained from the sea state program [Ref. 9]. Data input to 
the sea state program for Mariner Class Ship is shown in 
[Refs. 1, 10]. 

To see the ship's behaviour in sea state, the simulation 
program was run for ship speed 15 Knots, encounter angle 120 
degree, encounter frequency 0.64 radian/second, and sea 
state 8. First and second-order forces and moments were 
added into the surge, sway, and yaw equations that were used 
in the simulation program in Appendix C. Time responses of 
u, and 8 are presented in Figure 4.7, 4.8, and 4.9. In 
these simulations, optimal gain parameters found by the 
BOXPLX subroutine were used with autopilot. 

The first-order forces and moments are sinusoidal, 
therefore they create only oscillatory ship motions along 
the surge, sway, and yaw axes. The second-order forces and 
moments are constant, thus creating drift motions along 
these axes. 

It is apparent that as the encounter frequency 
increases, the effect of the first-order forces on ship 
motions decreases because of the high ship inertia ( acting 
as a low pass filter ) [Ref. 10, 21, 22]. An increasing 
encounter frequency means a decreasing wavelength- to- ship 
length ratio. The constant value of the second-order forces 
increases when the wavelength-to-ship length ratio 
decreases. Regarding the theory, forces and moments are 
significant for short wavelengths .[ Ref s . 9, 17] 
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Calm Water Turning Maneuver, Rudde: 
Yaw Rate vs. Time. 
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Water Turning Maneuver, 
Sway vs. Time. 
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Figure 4.2 
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Figure 4 . 4 



Calm Water Turning Maneuver, Rudder 
X vs. Y. 



25 deg. 
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Figure 4.5 
Y aw 



Calm Water Zig-zag Maneuver 
and Rudder vs. Time. 




Figure 4.6 



Calm Water Zig-zag Maneuver 
Speed vs. Time. 
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Figure 4. 7 
Encounter 



S . State = 8, Speed=15 Knots, D Type Cont . 
Angle=120 Deg., Encounter Frequency=0 . 64 
Yaw Degree vs. Time. 




Figure 4.8 
Encounter 



S. State = 8, Speed=15 Knots, D Type Cont. 
Angle=120 Deg., Encounter Frequency=0 . 64 
Rudder vs. Time. 
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Figure 4.9 S. State = 8, Speed=15 Knots, D Type Cont. 
Encounter Angle=120 Deg., Encounter Frequency=0 . 64 
Speed with Increasing RPM vs. Time. 
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V. COST FUNCTION 



It is accepted that the performance objective of the 
system in open- sea course keeping is " minimum added resis- 



fair-weather trial speeds for ships, the increase of resis- 
tance due to steering is often taken to be one percent of 
the towline resistance. Although this may offer little for 
improvement, when it is based on a fuel cost 0.5 percent of 
10 million dollars a year still equals 50000 dollars. 
Moreover, in moderate to bad weather the benefits of good 
steering rapidly become more evident. It is difficult to 
measure the changes of rate of fuel consumption due to 
different autopilots. A number of alternative cost functions 
have been suggested to describe the rate of fuel consumption 
[Refs. 23, 24, 25, 25, 27, 28]. 

When deriving a cost function for the autopilot, a 
requirement is that it must be convenient for ship board 
use. The cost function that has been commonly used in recent 
years is: 



0 

This is an approximate form of the exact cost function. 
A. EXACT FORM 

The performance criterion which characterizes the 
propulsion losses due to steering may be shown to be that 
derived from excess power consumption per unit distance 
caused by the added resistance due to steering. The added 
resistance due to steering can be related to the surge equa- 
tion where the total instantaneous surge relevant to 
steering is [Ref. 29]: 



tance due to steering " . For Routine predictions of the 



T 




( eqn 5.1) 
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( eqn 5.2) 



AX = (l/2)X w v 2 + (1/2)X S8 S 2 + (X vr +m)vr 
+ ( 1/2 ) ( X rr +mX G ) r 2 + X^vS 



This may be used to determine the energy losses related 
to steering in both calm water and waves. It is not conven- 
ient for ship board use, but may be used when finding fuel 
cost and comparing two controllers in fuel consumption. 

Then the exact form of the performance criteria is: 

T 

J = T lim ( 1/2T) f [(l/2)X vv v 2 + (1/2)X SS S 2 + (X vr +m)vr 

0 ~ 

+ ( 1/2 ) (X rr +mX G ) r 2 + X^vfl ] dt (eqn 5.3) 



The lower limit is written as zero, but while using a 
simulation program, the cost function will be calculated 
after tQ the transient response time of the system. 



B. APPROXIMATE FORM 

Accurate knowledge of the nonlinear coefficients in the 
exact form of the performance criterion is required for 
accurate results. In addition the criterion itself suffers 
from the disadvantage that the sway velocity is not measur- 
able in practice. 

How can the approximate form which is suitable for ship 
board be found (Refs. 10, 23, 24]? 

Since the sway velocity of the ship is small, the term 
including the square of the sway velocity can be neglected. 
It is also seen that X v ^v5 and ( 1/2 ) (X rr +mXQ) r 2 terms are 
small compared to others [Ref. 6]. After these assumptions. 
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the cost function for the Mariner Class Ship can be written 



as : 



= lim ( 1/2T ) 

!’-*• 03 



(l/2)X^^d + (X vr +m)vr ]dt ( egn 5.4) 



vr 



v and r are difficult to measure on ship board. vr can 
be defined in more convenient form for ship board use. 

Yawing and swaying of the vessel is assumed to stem from 
either self oscillations due to the steering system or the 
forced oscillations due to waves. The ship motions resulting 
from oscillatory seaway disturbance are oscillatory. 
Although these oscillations do not exactly follow a sinu- 
soidal pattern, they may be approximated as regular yawing 
of simple periodic form. Then yaw rate, sway velocity, and 
rudder angle can be represented as: 



r = r a Sin(cot+£ r ) 

v = v a Sin( a>t+<j) v ) ( eqn 5.5) 

<5 = 6' a Sin(wt+p£) 

Therefore : 



J = (1/2 )X SS S 2 Sin 2 ( W t+<^) 

+ (X vr +m) (v a r a /2) Cos(£ v -£ r ) (eqn 5.6) 

It is assumed that during low frequency oscillations 
Cos ( <£ v - <£ r ) - ■ 1 , i.e. , yaw rate and drift angle are in phase 
with one another. And writing jS=v/U as drift rate and 
assuming small amplitude oscillations around the pivot point 
p then, from Figure 5.1: 
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0 = (OP/R) = ( OP/L ) * ( L/R ) 



( eqn 5.7) 



where 



R = (U/r) 



Therefore : 



vr = ru 0 = ru( OP/L) (L/R) = r 2 OP 



( eqn 5.8) 



where 



r — \p — u)\p 



Then the expression for the cost function becomes: 



i 

J = lim ( 1/2T ) f [(1/2 )X SS 8 : 



+ ( X +m ) OP co 2 r// 2]dt (eqn 5.9) 



Then: 



J = lim ( 1/2T) 

T-*- co 



i 

/■ 



(\xp 2+ S 2 )dt 



( eqn 5.10) 



where 



(X vr + m) OPco 2 

X = 
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G/C 




Figure 5 . 1 



Geometry of 



Ship Turning. 
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T * T 
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T7 T TT7 r r, ^A T Q T !Vr rr ' T 

The accuracy of the ship's path is irrelevant on the 
open sea. The important thing is minimizing the propulsion 
power losses. These propulsion power losses can be caused by 
added drag due to steering and external forces and moments 
due to sea waves and wind acting on the ship. 

Because e and 5 values are easily measured, the 
optimum control parameters can be found by using the approx- 
imate form of the performance criterion. However, this does 
not provide evaluation of fuel consumption, nor does it 
permit comparison of different controller designs. 

Since ship operation may try to maintain constant power 
at the engine, or may attempt to maintain constant speed, 
autopilot design might be required for either condition of 
operation. 



M C * VT* IT V ^ T r T IT 

wWJJ — A. “VI » «L. 1 Itu 



POWER 



In calm water, the ship is allowed to proceed forward at 
a steady speed, and is free from lateral disturbance while 
undertaking part of a voyage. In this case the ship should 
take an elapsed time t, to cover a distance s, and speed u. 
In reality, the ship does encounter lateral disturbances 
which causes it to yaw at angular velocity r, sway sideways 
with velocity v, and rudder resistance is increased while 
keeping the course. This causes the ship to travel a longer 
path (As) than necessary, a longer time (At) is required 
than for the undisturbed case, and there is speed loss 
(Au) . 

The engine uses a certain amount of fuel in a given 
time : 



= W 



P 



net • 



m 



(eqn 6.1) 



Fuel Consumption = P . t 



(eqn 6.2) 



In the disturbed case: 



Fuel Consumption = P . t' 



(eqn 6.3) 



where 



t' = t+At 



The ship path can be simulated with the computer simula- 
tion program in Appendix D. And the time to travel the same 
distance along the x-axes can be determined for calm water 
and for the sea state case. 

In fact, to find how much fuel is consumed, the engine 
specifications must be known. But, when different control- 
lers are compared to determine the percent fuel savings 
in the same situations, there is no need for engine 
specifications . 

% Fuel Consumption P t^ “ P 

Ratio = x 100 (eqn 6.4) 



P*t '2 = Fuel consumption at t ' , elapsed time to travel 
the given distance along the x-axes while using controller 
#1( assume as an automatic optimal controller) 

P»t '2 = Fuel consumption at t^, elapsed time to travel 
the given distance along x-axes while using controller 
#2(assume as a PID controller) 




where 



then : 



% Fuel Consumption 
Ratio 




x 100 (eqn 6.5) 
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The sea state forces and moments, and added drag due to 
steering affect the' elapsed time to travel the given 
distance along the x-axes. This ratio gives the fuel 
consumption ratio including all sea state forces, moments, 
and added drag due to steering. But, the effect of added 
drag due to steering for different controllers is needed to 
compare fuel consumption. 

The other disadvantage is that it can be used only for 
small sea state disturbances, since the ship speed decreases 
too much in high sea state. It may be even sufficient to 
drive the ship backwards [Ref. 6). The vital effect of the 
decreasing speed is that this changes the values of the 
hydrodynamic coefficients of the nonlinear equations given 
in Table 2, 3, and 4 in Appendix B. The big speed differ- 
ences from steady-state speed, 15 Knots, makes the nonlinear 
coefficient values incorrect to represent the ship-steering 
dynamics . 

As a result; for small disturbances, i.e. , when the 
representation of ship- steering dynamics does not change, 
then this method may be used. But, for large sea distur- 
bances, the ship speed must not be allowed to decrease much 
to affect the representation of ship- steering dynamics in 
the simulation program. 



B. CONSTANT SPEED 

Using Table 6.1 [Ref. 6), it is seen that the added drag 
force due to increasing RPM is a linear function. The linear 
equation " added force vs. RPM " is found, and used in the 
simulation program to keep the speed constant. 

While attempting to maintain speed constant, it is more 
difficult to find the fuel consumption ratio, since P is 
varying. P cannot be factored out, as in Eqn. 6.5. 

How can a method to find the fuel consumption ratio 
without using engine specifications be found? 
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The added resistance due to steering is: 



AX = (l/2)X vv v 2 + (l/2)(X rr +mX G )r 2 + ( 1/2 ) X^ 2 



+ (X vr +m)vr + X v gv5 



(eqn 6.6) 



The total surge equation is: 



X total ~ x + x calm 



(eqn 6.7) 



where 




The fuel consumption ratio of added resistance due to 
steering to the total surge equation, assuming constant 
overall propulsive efficiency, [Refs. 2, 25] is: 



Using this ratio, it is possible to compare the fuel 
consumption of the Automatic Optimal Controller and with 
that of the PID. To find exact fuel consumption, engine 
specifications are still needed. This ratio will be used in 
Chapter 7 to compare the PID with the Automatic Optimal 
Controller . 



% Fuel Consumption 
Ratio 



Ax 



(eqn 6.8) 
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TABLE I 




SPEED AND ADDED FORCE VS 


. RPM 


RPM 


Speed 


Added Force 




( Knots ) 


(LB) 


55.5 


7.0 


13320 


56.8 


8.0 


16870 


58.1 


9.0 


21070 


59.6 


10.0 


25710 


61.2 


11.0 


31400 


62.8 


12.0 


37780 


64.6 


13.0 


45500 


66.5 


14.0 


53960 


68.6 


15.0 


63420 
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VII. COMPARISON OF PID AND AUTOMATIC OPTIMAL CONTROLLER 



The Automatic Optimal Controller offers the potential 
for minimization of propulsion losses due to steering in the 
open sea, and for removal of operator judgement from 
steering control of the ship. The optimal parameters which 
were found by using the BOXPLX subroutine in the Optimal 
Gain Program were shown in Table II [Ref. 1]. The approxi- 



mate form of the 
subroutine . 


cost 


function was used 


in 


the BOXPLX 


Time responses 


of \p 


and 5 are shown 


in 


figure 7 . 1 



through figure 7.10 for different encounter angles using the 
optimal parameter values in Table II . 

No provision for automatic adaptivity to either speed, 
load or seaway exists in the Universal Gyropilot, PID. Some 
adjustments of control parameters are possible, however, 
through an operator interface. Details of such adjustments, 
and the structure and parameter values of a PID controller 
are contained in Appendix A. 

Although a PID controller can come close to emulating 
this performance under specified internal and operator 
settings, the operator has no means to judge steering 
performance other than by observation of the course 
recorder. The normal tendency of the operator is to change 
the external controls available to him to bring about the 
reduction in heading error. It greatly increases the 
propulsion losses due to steering. 

The goal in this chapter is to compare an Automatic 
Optimal Controller with a well tuned PID controller in fuel 
savings. Now the question is how the Well -Tuned PID can be 
defined. There are two ways to answer this question. One way 
is to define a well tuned PID as an Optimal PID Controller 
and to use the same performance criterion used in finding 
optimal parameters of a D Type Automatic Controller, i.e.. 
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the minimization of added resistance due to steering. The 
optimal control parameters of a PID in different sea state, 
speed, encounter angle, and encounter frequencies were found 
by using the simulation program to compute the cost function 
of all combinations of internal and external control parame- 
ters, making the modifications in the simulation program 
shown in Appendix C. The results were shown in Table III . 

Time responses of \p and 5 are shown . in Figure 7.11 
through 7.20 for different encounter angles using the 
optimal parameter values in Table III . As it is seen from 
the figures, the Automatic Optimal Controller and PID with 
optimal parameters after the transient time makes the ship 
yaw sinusoidally around the commanded course. To see their 
effects in added resistance due to steering, and in fuel 

i 

savings, the comparison of them was shown in Table IV . The 
added resistance due to steering and fuel consumption ratio 
are almost the same, since the same performance criterion 
was used, and all internal and external controls of PID are 
adjusted, not only the external controls. 

The PID with optimum parameters cannot be used as a well 
tuned PID, since the operator could not find the optimum 
parameters by himself to minimize the added resistance due 
to steering and also could not change the internal control 
parameters. But it is clear that if the optimum control 
parameters are found in the same way as for the D Type 
Controller, it gives the same fuel savings as an optimal D 
Type controller. 

The other way is to approach the question as an oper- 
ator. There are no means to adjust the external controls 
available to the operator to minimize the propulsion losses 
due to steering other than by observation of the course 
recorder. The operator can adjust the external controls to 
reduce the heading error, the only system output available 
to him. The external controls available to the operator are 
Weather Adjust Gain, Rudder Multiplier Gain, and Rudder Rate 
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Multiplier. To reduce the heading error, the Weather Adjust 
Gain is set to "0", and to get the fast response of rudder, 
the Rudder Multiplier Gain is at "3", and the Rudder rate 
Multiplier setting is "l". These correspond to K-^ is "1", K 2 
is "3", and the value of T-^ is half of the internal control 
setting . 

The internal controls T^ T 2 and T 3 are not known for 
the Mariner Class Type ship. These internal control parame- 
ters are found by running the simulation program used in 
finding optimum parameters for the PID with and K 2 set to 
"1", and the proportional part of the PID is omitted in the 
control law of PID, Egn. A.l. The results are following: 

= 25 
T 2 = 15 

t 3 = 1000 

Time responses of \p and § are shown in Figure 7.21 
through 7.30. As is seen in Figure 7.22, 7.24, 7.26, 7.28, 
and 7.30, the well tuned PID uses a larger rudder angle to 
reduce the heading error quickly. The rudder angle also 
oscillates in large magnitudes. It tries to get the ship to 
the command heading angle as fast as possible. It increases 
the cost. So consequently the fuel cost and added resis- 
tance due to steering increase. The comparison in fuel 
savings and in added resistance due to steering is shown in 
Table V . Although the optimal internal control parameters 
are used in a well tuned PID, there are differences in fuel 
consumption ratio and in added resistance due to steering, 
especially at quartering sea waves, i.e., 030-060 and 
120-150 degrees. The Automatic Optimal Controller provides 
fuel savings in excess of 0.5 % over a well tuned PID when 
operating at the design speed at random headings in sea 
states . 
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The internal control parameters of a well tuned PID in 
the deterministic model are represented by optimal values. 
They might not be the actual values used for a Mariner Class 
ship. The differences in fuel consumption ratio and in added 
resistance due to steering should increase if actual values 
are not equal to the optimal internal control parameter 
values used in the simulation computer program. It is also 
possible that the different operators can set the Rudder 
Multiplier Gain and Rudder Rate Multiplier in different 
ways. The fuel consumption ratio and added resistance due to 
steering also changes in relation to these values. 
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TABLE II 

OPTIMAL PARAMETERS OF D TYPE CONTROLLER 

Sea State = 8 

Ship Speed = 15 Knots 

Wave Frequency = 0.53 

Encounter 

Angle K T z T p K d 



030 


0.69 


0.01 


34.38 


569.14 


060 


6.70 


0.01 


73.47 


23 . 16 


090 


7.45 


0 . 19 


60.58 


990.44 


120 


2.77 


1.75 


35.79 


29.76 


150 


2 . 41 


0.01 


6.29 


12 . 42 



TABLE III 

OPTIMAL PARAMETERS OF PID CONTROLLER 



Sea State = 8 
Ship Speed = 15 
Wave Frequency : 


KNOTS 
= 0.53 








Encounter 

Angle 


K 1 


k 2 


T 1 


T 2 


T 3 


030 


0.33 


1 . 5 


10. 


15 . 


200 


060 


1 . 


1.5 


25 . 


10. 


200 


090 


1 . 


1 . 


25. 


15. 


200 


120 


1 . 


1 . 


2.5 


15. 


200 


150 


1 . 


2. 


2.5 


15. 


200 
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TABLE IV 



COMPARISON OF OPTIMUM PID WITH AUTOMATIC OPTIMAL 

CONTROLLER 



Sea State = 8 
Ship Speed = 15 Knots 
Wave Frequency = 0.53 
Added Resistance (libre) 

Encounter Angle 030 deg. 





D Type 


Optimum PID 


Difference 


Added Resistance 
Due to Steering 


5547 


5892 


345 


Fuel Consumption 
Ratio ( % ) 


6. 106 


6.552 


0.446 



Encounter Angle 060 deg. 





D Type 


Optimum PID 


Difference 


Added Resistance 
Due to Steering 


52422 


52734 


312 


Fuel Consumption 
Ratio ( % ) 


33 . 947 


34.080 


0.133 



Encounter Angle 090 deg. 





D Type 


Optimum PID 


Difference 


Added Resistance 
Due to Steering 


108930 


109250 


320 


Fuel Consumption 
Ratio ( % ) 


59.749 


59.928 


0.179 


Encounter Angle 120 deg. 






D Type 


Optimum PID 


Difference 


Added Resistance 
Due to Steering 


51377 


51857 


210 


Fuel Consumption 
Ratio ( % ) 


28. 179 


28.442 


0.179 


Encounter Angle 150 deg. 






D Type 


Optimum PID 


Difference 


Added Resistance 
Due to Steering 


14756 


15221 


465 


Fuel Consumption 
Ratio ( % ) 


8.052 


8.307 


0.255 
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TABLE V 



COMPARISON OF WELL TUNED PID WITH AUTOMATIC OPTIMAL 

CONTROLLER 



Sea State = 8 

Ship Speed = 15 Knots 

Wave Frequency = 0.53 



Encounter Angle 030 deg. 





D Type 
Cont . 


Well Tuned 
PID 


Difference 


Added Resistance 
Due to Steering 


5547 


7178 


1631 


Fuel Consumption 
Ratio ( % ) 


6.106 


7.761 


1 . 655 



Encounter Angle 060 deg. 





D Type 
Cont. 


Well Tuned 






, P ID 


Difference 


Added Resistance 
Due to Steering 


52422 


53894 


1108 


Fuel Consumption 
Ratio ( % ) 


33.947 


34.340 


0.393 



Encounter Angle 090 deg. 





D Type 
Cont . 


Well Tuned 






PID 


Difference 


Added Resistance 
Due to Steering 


108930 


109267 


337 


Fuel Consumption 
Ratio ( % ) 


59 . 749 


59.934 


0.185 



Encounter Angle 120 deg. 





D cSr 


Well Tuned 






PID 


Difference 


Added Resistance 
Due to Steering 


51377 


53561 


2184 


Fuel Consumption 
Ratio ( % ) 


28.179 


29.375 


1 . 196 



Encounter Angle 150 deg. 





D Type 
Cont. 


Well Tuned 






PID 


Difference 


Added Resistance 
Due to Steering 


14756 


15747 


991 


Fuel Consumption 
Ratio ( % ) 


8.052 


8.593 


0.541 
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Figure 7.1 S. State=8, D Type Cent., Speed=15 Knots 
Encounter Angie=C3Q, Encounter Frequency=C . 33 
Yaw Degree vs. Trrr.e. 




Figure 7.2 S. State=8, D Type Cont., Speed=15 Knots 
Encounter Angle=030, Encounter Frequency=0 . 33 
Rudder vs. Time. 
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Figure 7.4 S. State=8, 
Encounter Angle=060, 



D Type Cont . , Speed=15 Knots 
Encounter Frequency=0 . 43 



Rudder vs. Time. 




Figure 7.5 5. State=8, D 

Encounter Angie=090, E 



ype Cont. , Speea=15 Knots 
cutter Frequency=0 . 53 



Yaw Degree vs. Time. 




Figure 7.6 S. State=8, 
Encounter Angle=090, 
Rudder 



D Type Cont., Speed=15 Knots 
Encounter Frequency=0 . 53 
vs. Time. 
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Figure 7.8 S. State=8, 
Encounter Angle=120, 



D Type Cont . , Speed=15 Knots 
Encounter Frequency=0 . 64 



Rudder vs. Time 



in 








1 

' L 


i! 

kj 

0 i 11 ' 






IT > - 
















n ri, 






CJ)r\j *1 
Q) 

~D 


in 




Siib'jS 


2 °" 




ii,J 




az 

>-< 

in 


y 1 

i* 


H 

SEL3:i 




r\3 ~ 
» 

ip 


ii 


W 





0 . CO 



h J L* Jl. «... J . * >1 * ii L, i t t ■ « .. i it 1 ■. j / l(il ^ ) 1 



a 






iHirlWljiiii; 






ftp 



Mliii iMlMiul 






j/ * U.J- . Uj 

TIME isec. 



’igure 7.9 S. State=3, D Type Con t., Speed=15 Kr.c 
Encounter Angie=150, Encounter F requer.cy-0 . 5" 
Yaw Dearee vs. Time. 




Figure 7.10 S. State=8, D Type Cont. , Speed=15 Knots 
Encounter Angle=150, Encounter Frequency=0 . 67 
Rudder vs. Time. 
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Figure 7 . 12 

Encounter 



S. State=8, Optimal FID, Speed=15 Knots 
Angle=030 , Encounter Frequency=0 . 33 
Rudder vs. Time. 
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Figure 7.14 
Encounter 



S. State=8, Optimal PID, Speed=15 Knots 
Angle=060, Encounter Freguency=0 . 43 
Rudder vs. Time. 
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Figure 7.16 

Encounter 



S. State=8 / Optimal FID, Speed=15 Knots 
Angle=090 , Encounter Frequency=0 . 53 
Rudder vs. Time. 
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Figure 7. 
Zr.cou 



S. State—S, Optimal FID, Speed=15 Kn 
Angie=12C, Encounter Frequency=D . 54 
Ya\; Decree vs. Time. 




Figure 7.18 

Encounter 



S. State=8 / Optimal PID, Speea=15 Knots 
Angle=120, Encounter Freguency=0 . 64 
Rudder vs. Time. 
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Figure 7 . 19 



Encou 



e r 



S. State=8, Optimal FID, Speed=15 Knots 
Ancie=l 5 0 , Encounter Erecuencv=0 . 57 
Yaw Degree vs. Tine. 




Figure 7.20 

Encounter 



S. State=8, Optimal FID, Speed=15 Knots 
Angle=150, Encounter Frequency=0 . 67 
Rudder vs. Time. 
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Figure 7.21 S. State=8, Well Tuned FID, Speed=15 Knots 
Encounter Angie=C30, Encounter Frequencv=Q . 33 
Yaw Degree vs. Time. 




Figure 7.22 S. State=8, Well Tuned PID, Speed=15 Knots 
Encounter Angle=030, Encounter Frequency=0 . 33 
Rudder vs. Time. 
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Figure ".23 S. State=S, Well Tuned FID, Speed=15 Knots 
Encounter Angie=Q60, Encounter Frequency=0 . 43 
Yaw Degree vs. Tire . 




Figure 7.24 S. State=8, Well Tuned PID, Speed=15 Knots 
Encounter Angle=060, Encounter Frequency=0 . 43 
Rudder vs. Time. 
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Figure 7.26 S. State=8, Well Tuned PID, Speed=15 Knots 
Encounter Angle=090, Encounter Frequency=0 . 53 
Rudder vs. Time. 
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Figure 7.28 S. State=8 / Well Tuned 
Encounter Angle=120, Encounter 
Rudder vs. Time. 
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Figure 7.29 S. State=2, Weil Tuned FID, Speec=15 Knees 
Encounter Angle=150, Encounter F requeue y=0 . 67 
Yaw Degree vs. Tine. 





Figure 7.30 S. State=8, Well Tuned PID, Speed=15 Knots 
Encounter Angle=150, Encounter Frequency=0 . 67 
Rudder vs. Time. 
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VIII. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

An A Type controller was used in finding optimal parame- 
ters to minimize the propulsion losses due to steering in 
previous theses [Refs. 8, 10, 21, 22]. But to prevent the 
drifting effect of second-order forces and moments on the 
ship's motions, a D Type controller was used in designing 
the Automatic Optimal Controller. The only difference 
between them is the integrator part omitted in m the A Type. 

Savings in fuel for an Automatic Optimal Controller, an 
Optimal PID, and a Well-Tuned PID were compared by using a 
new approach explained in Chapter 6 that does not use engine 
specifications. It can be used in finding any ship's fuel 
savings by finding the added resistance due to steering. 

The optimal internal and external settings of a PID to 
minimize the propulsion losses due to steering are found by 
the modified simulation program shown in Appendix C, since 
there is not enough information about the settings .of 
control parameters of PID for the Mariner Class ship. The 
same performance criterion explained in Chapter 5 was used 
in finding the optimal parameters for both of them. As a 
result, it is found that the optimal PID is as - good as 
Automatic Optimal Controller. They provide almost the same 
fuel savings in the same situations. 

The well tuned PID was defined as a PID controller to 
reduce the heading error as fast as possible and to get fast 
response of the rudder to do this. Only external control 
settings (explained in Appendix A) are available to the 
operator. The optimal internal control settings were used in 
simulations . The Automatic Optimal Controller provides fuel 
savings in excess of 0.5 % over a well tuned PID when oper- 
ating at sea state, 8, speed, 15 Knots, and wave frequency, 
0.53 radian/seconds. 
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B. RECOMMENDATIONS FOR FUTURE STUDIES 

Additional work should be done to compare the Adaptive 
Automatic Optimal Controller with a well tuned PID for 
different sea states, speeds, and encounter frequencies. 

In this thesis the well tuned PID was used as explained 
in Chapter 7. If more information is available for the PID, 
especially if a precise definition can be found for the well 
tuned PID, then the comparison in fuel savings with these 
well tuned PID settings will provide better insight into the 
value of the Optimal Controller. 

In this thesis, surge, sway, and yaw equations of motion 
were used in simulations. The roll equation of motion may 
be added to the simulation program and the optimal parame- 
ters of the Automatic Optimal Controller may be found by 
using the BOXPLX minimization subroutine. Then the compar- 
ison in fuel savings will be more realistic. 
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APPENDIX A 



PID CONTROLLER 

The standard form of a PID controller is based on 
proportional-plus-integral -plus-derivative control. No 
provision for automatic adaptivity of the controller to 
speed or seaway effects exists in the PID. It relies on 
adjustment of the control parameters by the operator. 

1. GENERAL DESCRIPTION 

An autopilot commonly found aboard merchant ships has 
the characteristics of a PID controller. The PID control law 
is described by [Refs. 23, 24]: 



Tis 



= K-| [ Ko(l + 



(1+T 2 s )' 



) + 



T 3 s 



] \p i 



where the gains are determined by operator interaction in 
the following manner: 

1 . Weather Adjust Gain 

= 1/3 within "weather adjust" zone 
= 1 outside "weather adjust" zone 
The weather adjust zone is variable from 0.3 to 5.0 
degree heading error. Panel is marked "0" to "5". 

2 . Rudder Multiplier Gain 

K 2 = 1 to 3 , continuously adjustable (panel 

control ) 
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3 . Internal Control Settings 

The time constants T-^, T 2 , and T 3 are preset for the 
specific ship, with the exception that may be reduced by a 
factor of 2 by operator interaction. The ranges of possible 
gyropilot time constants are: 

1. Derivative Time Constant 

T^ = 2.5, 5, 10, 25, 50, 100, 150, 200 second 
"Rate Mult" control reduces the time constant by a 
factor of 2 when at minimum setting. 

2. Derivative Filter Time Constant 

T 2 = 1.5, 2.5, 5, 10, 15 second 

3 . Integral Time Constant 

T 3 = 200, 400, 600, 1000 second 

2. APPLIED TO SIMULATION PROGRAM 

Figure A.l corresponds to controller PID which has the 
form [ ref . 22 ] : 




1 



K x [ K 2 (l + 



This can be written in simulation program as this: 



DX2 = ( YAWE - X2 ) / T 2 
DX3 = ( X2 - X3 ) / T 2 

D = ( K 1 *K 2 *T 1 *DX3 ) + ( K 1 *K 2 *YAWE) + X4 



Where 



X2 = X2 + ( DX2 * DELT ) 

X3 = X3 + ( DX3 * DELT ) 

X4 = X4 + ( YAWE * (K 1 /T 3 ) ) * DELT 
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Figure A.l PID Controller in Simulation Program. 
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APPENDIX B 



NONLINEAR EQUATIONS OF MOTION AND COEFFICIENTS TABLES 
Nonlinear Equations of Motions: 



X-Eq. : 


(m - X^)U = 


f 1 (u / v, r , 5 ) 






* 


9 


f 2 (u,v, r,<5) (eqn B. 1) 


Y-Eq. : 


(m - Y^)v + (mX G - Y^)r = 


N-Eq. : 


(mX G -N^)v + (I 2 - N£ ) r = 


f 3 (u, v, r,5) 



Where 

f 1 (u,v,r < 5)= X* + X* u Au + (l/2)X uu Au 2 + ( 1/6 ) X uuu Au 3 + 
(l/2)X vv v 2 + ( (l/2)X rr +mX G )r 2 + ( 1/2 )Xo 5 «$ 2 + 
(1/2 )X vvu v 2 Au + (l/2)X rru r 2 Au + ( 1/2 )X 55u ,5 2 Au + 
(X vr +m)vr + X V(5 v< 5 + X^r <5 + X vru vrAu + 
x V(5u v 5Au + x r5u r,5Au 
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f '2 ( u ' v ' r, 5 )= Y* + Y* u Au + Y* 1JU Au 2 + Y v v + (l/6)Y vvv v 3 + 
(l/2)Y vrr vr 2 + (1/2 )Y v55 v§ 2 + Y vu vAu + 
(1/2 )Y vuu vAu 2 + (Y r -mu)r + (l/6)Y rrr r 3 + 
(l/2)Y rvv rv 2 + (l/2)Y r(S(5 r(5 2 + Y ru rAu + 
(l/2)Y ruu rAu 2 + Y S S + (1/6)Y SSS S 3 
d/2)Y Svv Sv 2 + y 5u 5Au + (1/2)y 5uu 5Au 2 + 
Y vr5 vr<5 + (l/2)Y 5rr 5r 2 

f 3 ( u / v, r ,S)= N* + N* u Au + N* uu Au 2 + N v v + (l/ 6 )N vvv v 3 + 
(l/2)N vrr vr 2 + (1/2 )N v5(5 v 5 2 + N^vAu + 
(l/2)N vuu vAu 2 + (N r -mX G u)r + (l/6)N rrr r 3 + 
(l/2)N rvv rv 2 + (l/2)N r55 r5 2 + N ru rAu + 
(l/2)N ruu rAu 2 + N 5 S + (1/6)N 555 5 3 + 

( 1/2 )Ng vv <5v 2 + N 5 u 5Au + d/2)N fiuu 5Au 2 + 
N V rfiVr5 + ( 1/2 )N 5rr 5r 2 
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TABLE VI 




COEFFICIENTS OF 


X- EQUATION 




Taylor 
Expansion 
and Dynamic 
response 
Terms 


Identifier 

in 

Fortran 

Program 


Nondim. 

Factor 


Nondim. 
Coeff . 

*10 5 


(m-X^) 


Cl 


( 1/2 ) p LBP 3 


840.0 




XI 


( 1/2 )p LBP 2 S- 


-120.0 


d/2)^ 


X2 


( 1/2 ) p LBP 2 


45 . 0 


d/ 6 )Xuuu 


X3 


( 1/2 )p LBP 2 /S 


-10.3 


( 1/2 )X VV 


X4 


( 1/2 ) p LBP 2 


-898.8 


(l/2)X rr +mX G 


X5 


(1/2 )p LBP 4 


18.0 


(V6)X 5(5 


X6 


( 1/2 ) p LBP 2 S 2 


-94-8 


( 1/2 ) X wu 




( 1/2 )p LBP 2 /S 




<l/2)X rru 




( 1/2 ) p LBP 4 /S 




(V2)X 55u 




( 1/2 )p LBP 2 S 




(X vr +m) 


X7 


( 1/2 ) p LBP 3 


•798.0 


x v5 


X8 


( 1/2 ) p LBP 2 S 


83.2 


X r5 




( 1/2 )p LBP 3 S 


0.0 


y 

vru 




( 1/2 ) p LBP 3 /S 








( 1/2 )p LBP 2 




X r6u 




( 1/2 ) p LBP 3 




X* 


XO 


( 1/2 )p LBP 2 S 2 


0.0 
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TABLE VII 




COEFFICIENTS OF 


Y- EQUATION 




Taylor 
Expansion 
and Dynamic 
response 
Terms 


Identifier 

in 

Fortran 

Program 


Nondim. 

Factor 


Nondim. 
Coef f . 

* 10 5 


(m-Y$) 


C2 


( 1/2 ) p LBP 3 


1546.0 


(mX G -Y*) 


C3 


( 1/2 )p LBP 4 


-8.6 


Y v 


Y1 


( 1/2 )p LBP 2 S 


-1160.4 


< 1 /6)Y VVV 


Y2 


(1/2 )p LBP 2 /S 


-8078.2 


(l/2)Y vrr 




(1/2 )p LBP 4 /S 


0.0 


( 1 /2)Y v55 


Y3 


( 1/2 )p LBP 2 S 


-3.8 


( Y r -mu) 


Y4 


(1/2 )p LBP 3 S 


-499.0 


(l/6)Y rrr 




(1/2 )p LBP 5 /S 


0.0 


(l/2)Y rvv 


Y5 


( 1/2 ) p LBP 3 /S 


15356.0 


( !/ 2 ) Y r§5 




( 1/2 )p LBP 3 S 


0.0 


y 5 


Y6 


( 1/2 )p LBP 2 S 2 


277.9 


( 1 /6)Y SSS 


Y7 


( 1/2 ) p LBP 2 S 2 


-90.0 


(V2)Y 5vv 


Y8 


(1/2 )p LBP 2 


1189.6 


d/2)Y drr 




( 1/2 ) p LBP 4 


0.0 


y 5u 




( 1/2 ) p LBP 2 S 


0.0 


(V2)Y 5uu 




(1/2 )p LBP 2 




Y vr5 




( 1/2 )p LBP 3 


0.0 


Y* 


YO 


( 1/2 )p LBP 2 S 2 


-3 . 6 


Y* 

U 




(l/2)p LBP 2 S 


0.0 
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TABLE VIII 




COEFFICIENTS OF 


N- EQUATION 




Taylor 
Expansion 
and Dynamic 
response 
Terms 


Identifier 

in 

Fortran 

Program 


Nondim. 

Factor 


Nondim . 
Coef f . 

*10 5 


(mX G -N«) 


C5 


( 1/2 )p LBP 4 


-22 . 7 


(V N r) 


C4 


( 1/2 )p LBP 5 


82.9 


N v 


N1 


( 1/2 )p LBP 3 S 


-263.5 


(1/6)N WV 


N2 


( 1/2 )p LBP 3 /S 


1636.1 


(l/2)N vrr 




( 1/2 ) p LBP 5 /S 


0.0 


(1/2)N v65 


N3 


( 1/2 )p LBP 3 S 


12.5 


(N r -mX G u) 


N4 


( 1/2 )p LBP 4 S 


-166.0 


(l/6)N rrr 




( 1/2 ) p LBP 6 /S 


0.0 


(l/2)N rvv 


N5 


( 1/2 )p LBP 4 /S 


-5483.0 


(V 2 )^ r88 




( 1/2 ) p LBP 4 S 


0.0 


n 6 


N6 


( 1/2 ) p LBP 3 S 2 


-138.8 


d/6)N 5(55 


N7 


( 1/2 )p LBP 3 S 2 


45 . 0 


(V2)N 6w 


N8 


( 1/2 )p LBP 3 


-489.0 


(V2)N 5rr 




(l/2)p LBP 4 S 


0.0 


n 5u 




( 1/2 )p LBP 3 S 


0.0 


d/2)N vr5 




( 1/2 )p LBP 4 


0.0 


N* 


NO 


( 1/2 )p LBP 3 S 2 


2.8 


N* 

N u 




( 1/2 )p LBP 3 S 


0.0 
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APPENDIX C 



SIMULATION PROGRAM 

C*** DECLERATIONS OF VARIABLES ************************** 



REAL 


MH( 10000) ,T( 10000) 


REAL*8 


L,L2,L3,L4,L5,L6 


REAL* 8 


X, XDOT, Y, YDOT, U, UDOT, V, VDOT, YAW, R, RDOT 


REAL*8 


XI , X2 , X3 , X4 , X5 , X6 , X7 , X8 


REAL*8 


Y0,Y1,Y2,Y3,Y4,Y5,Y6,Y7,Y8 


REAL*8 


NO , N1 , N2 , N3 , N4 , N5 , N6 , N7 , N8 , TEM, TEMP , TEM1 


C 


, TEMPI 


REAL*8 


Cl , C2 , C3 , C4, C5 , FI , F2 , F3 , F2X, F2Y, M2Z, FXX , FYY 


C 


, MZZ 


REAL*8 


RO, DELT , S , DU, U1,K,Z1,Z2,P1,P2, KD , RDI F 


REAL*8 


DYAWE , YAWE , YAWC , I SR , I SE , TD I FF , LAMDA , RPM , RES 


C 


, EHP 


REAL*8 


SI , S2 , DS1 , DS2 , D, LOST , YAWDEG, DDEG, S3 


REAL*8 


MASS , I Z , XG , YVDOT , NVDOT , YR , YRDOT , NR , NRDOT , FX 


C 


, FY, MZ 


REAL*8 


RX, RY, RZ , TX, TY , TZ , WA, WE , RXR, RYR, RXI , RYI , MZR 


C 


, MZI 


REAL*8 


MAGMAX, MAGMIN , TIME , ETIME , DELST1 , DELST2 , DELSTE 


REAL*8 


BARFUE , LBFUEL 


REAL*8 


TOTFR1 , TOTFR2 , PERC1 , PERC2 



C* * *WHEN PID CONTROLLER IS USED , ADD THE FOLLOWINGS ****** 
0*** VARIABLES ****** 

C REAL*8 PIDV2 , PIDV3 , PIDV4, PIDDV2 , PIDDV3 , PI DTI , PIDT2 

C C , pIDT3 , PIDK1 , pidk2 

C***WHEN OPTIMAL PID PARAMETERS ARE WANTED TO FIND,***** 
C *** ADD THE FOLLOWING VARIABLES AND GIVEN DATA VALUES*** 

C REAL*8 KK1 ( 2 ) , KK2 ( 5 ) , TT1 ( 9 ) , TT2 ( 5 ) , TT3 ( 4 ) 

C DATA KK1( l)/0. 33/1 ./ 

C DATA KK2 ( 1 ) , KK2 ( 2 ) , KK2 ( 3 ) , KK2 ( 4 ) , KK2 ( 5 ) /I . , 1 . 5 , 2 . 
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C C ,2.5,3./ 

C DATA TT1 ( 1 ) , TT1 '( 2 ) , TT1 ( 3 ) , TT1 ( 4 ) , TT1 ( 5 ) , TT1 ( 6 ) 

C C , TT1 ( 7 ) , TT1 ( 8 ) , 

C C TTl(9)/2 .5,5 . , 10. , 15. ,25. , 50. , 100. , 150. ,200./ 

C DATA TT2 ( 1 ) , TT2 ( 2 ) , TT2 ( 3 ) , TT2 ( 4 ) , TT2 ( 5 ) /I . 5 , 2 . 5 , 5 . 

C C ,10. ,15./ 

C DATA TT3( 1) , TT3 (2 ) ,TT3 (3 ) , TT3 (4)/200. , 400. , 600. 

C C .,1000./ 

c********************************************************** 

c 

c 

C*** THESE disspla statements are used only when ****** 

c*** RUNNING TO GET GRAPHS. THEY ARE NOT USED ****** 

C*** WHEN FINDING OPTIMAL PARAMETERS FOR PID. ****** 

CALL TEK618 

CALL PAGE (11. 0,8. 5) 

CALL NOBRDR 
CALL HEIGHT (0.' 15 ) 

CALL MX1ALF ( 'STANDARD ' , ' ) 

CALL MX2 ALF ( ' L/CSTD 

Q************************************************** 

C 

c*** USE THIS PART TO FIND OPTIMAL PID PARAMETERS ****** 
C DO 11 IK1=1 , 2 

C DO 22 IK2=1 , 5 

C DO 33 IT1=1 , 9 

C DO 44 IT2=1 , 5 

C DO 55 IT3=1 , 4 

C 
C 

C*** WHEN FINDING OPTIMAL PID PARAMETERS, ****** 

C*** USED ETIME=1000 . ****** 

ETIME=2000. 

TIME=0 . 

ITIME=0 
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IC0UNT=1 



C 

C 

C* ** THESE ARE OPTIMAL PARAMETERS OF AUTOMATIC OPTIMAL *** 
C*** CONTROLLER FOUND BY RUNNING OPTIMAL GAIN PROGRAM *** 
C*** IN REFERANS 7 , 9 , 21 , 22 , 23. *** 

K=0. 6983 
Z1=0 . 0102 . 

Pl=34 . 3839 
KD=569 . 1469 

c*** here, you can use well tuned pid parameters or **** 

C*** OPTIMAL PID PARAMETERS FOUND BY RUNNING THIS **** 

c *** SIMULATION PROGRAM USING THE PART FOR FINDING **** 

C*** OPTIMAL PID VALUES. **** 

C P IDK1=1 . 0 

C P IDK2=3 . 0 

C PIDT1=25. 

C PIDT2— 15 . 

C P IDT3=200 . 

C 

C*** INITIALIZE ALL THE NECESSARY VARIABLES *************** 
C INITIALIZE THE COST FUNCTION 
I SE=0 . 0 
I SR=0 . 0 
TD I FF=0 . 0 
LAMDA=5 . 48 

C X , XDOT , Y , YDOT ARE FIX COORDINATES ON EARTH 

X=0.0 
Y=0 . 0 
XDOT=0 . 0 
YDOT=0 . 0 

C U , UDOT , V , VDOT ARE FIX COORDINATES ON SHIP 

V=0.0 
UDOT=0 . 0 
VDOT=0 . 0 
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F1=0 . 0 
F2=0 . 0 
F3=0 . 0 
YAW=0 . 0 
YAWE=0 . 0 
R=0 . 0 
RDOT=0 . 0 

RPM=68 .6 . 

RES=0 . 

EHP=0 . 0 
BARFUE =0.0 
LBFUEL =0.0 
RO=l . 9876* . 5 
G=32 . 174 
YAWE=0 . 0 
TEM=0 . 0 
TEMP=0 . 0 
TEM1=0. 0 
TEMP 1=0. 0 
DS1=0 . 0 
DS2=0 . 0 
S1=0. 0 
S2=0 . 0 
P IDDV2=0 . 0 
P IDDV3=0 . 0 
P IDV2=0 . 0 
PIDV3=0 . 0 
P IDV4=0 . 0 
SMALL=99999999 . 9 

C D = RUDDER ANGLE 
D=0. 0/57. 296 

C ORDERED SPEED IN FEET/SEC 

C 15. *1.689 FT/SEC=15 KNOTS 
Ul=15. *1.689 

C AT STEADY STATE ACTUAL SPEED (U) = COMMAND SPEED (UC) 
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U=U1 
L=528 . 

L2=L**2 

L3=L*L*L 

L4=L*L3 

L5=L*L4 

L6=L*L5 

C FORCES IN X, Y DIRECTION COMPUTED IN FORCES 

C MOMENTS IN Z 

FX=0 . 

FY=0 . 

MZ=0 . 

C 

C 

C*** SEA STATE FORCES AND MOMENTS FOUND BY RUNNING***** 

THE JIM CASS'S SEA STATE PROGRAM IN REFERANS 8. ***** 

C*** FIRST ORDER FORCES AND MOMENT 
RXR= - . 12887D4 
RXI=- . 90512D5 
RYR= . 74476D3 
RYI= .11220D6 
MZR=- . 20222D8 
MZI=- . 10028D7 

C*** SECOND ORDER FORCES AND MOMENT 
F2X= . 90755D-3 
F2Y= . 16626D-4 
M2Z= . 96895D-4 

C 

C*** FINDING FORCES AND MOMENT MAGNITUDES AND PHASES ***** 
RX=(RXR**2+RXI**2 )** . 5 
RY= ( RYR**2 + RYI **2 ) ** . 5 
RZ= ( MZR* * 2 +MZ 1**2)** . 5 
TX=DATAN2 ( RXI , RXR ) 

TY=DATAN2 ( RYI , RYR ) 

TZ=DATAN2 ( MZ I , MZR ) 
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FXX=2 . *RO*G*L*WA* *2*F2X 
FYY=2 . *R0*G*L*WA**2*F2Y 
MZZ=2 . *RO*G* ( L**2 ) * ( WA**2 ) *M2Z 
C 
C 

C*** SEA STATE AND SHIP SPECIFICATIONS ******************* 
C SIGNIFICANT WAVE HEIGHT: SEA STATE 8-19 FEET, 7-15 

C . 6-11, 5-7 

WA=19 . 

C ENCOUNTER FREQUENCY :( WHEN ENCOUNTER ANGLE IS 30) 

WE=0 . 33 

C ADDED MASS AND ADDED INERTIA TERMS: 

MASS=. 11155D+07 
IZ=0 . 000392*RO*L5 
XG=- 11.96 
C 
C 

200 CONTINUE 
C INPUT YAW COMMAND 

YAWC= 0.0/57.296 
C 

C TO PREVENT THE SHIP TURNING MORE THAN ONE DEGREE 
C PER SECOND AT THE BEGINNING OF THE PROGRAM SINCE 
C THE SHIP' CANNOT BE IN HIGHER SEA STATE INSTANTLY. 

C THERE MUST BE TIME PERIOD C FROM GOING SEA STATE 
C ONE TO SEA STATE EIGHT. 

TEM1=YAW- TEMPI 

IF (TEM1.GT. 1.0/57.296) YAW=TEMP1+1 . 0/57 . 296 
IF(TEM1 . LT. -1.0/57.296) YAW=TEMP1 - 1 . 0/5.7 . 296 
TEMP1=YAW 
C 

C ERROR SIGNAL TO DRIVE RUDDER (YAW ACTUAL - YAW ORDERED) 
YAWE=YAW - YAWC 
C 

S=( (U*U)+(V*V) ) ** . 5 
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DU=U-U1 

c * * * I p D type AUTOMATIC OPTIMAL CONTROLLER IS USED,***** 
C***USE THESE ***** 

DS1= ( YAWE - SI) /PI 
D=( SI + DS1*Z1)*K+S2 
C 

c ***if OPTIMAL OR WELL TUNED PID IS USED TO GET GRAPHS*** 
C*** USE THIS PART. *** 

C P I DDV2 = ( YAWE - PIDV2 ) / PIDT2 

C PIDDV3= (PIDV2- PIDV3 ) / PIDT2 

C D = (PIDK1*PIDK2*PIDT1*PIDDV3 ) + ( P IDK1*PIDK2 *YAWE ) 

C C + PIDV4 
C 

C***WHEN FINDING OPTIMAL PID PARAMETERS, USE THIS PART.**** 
C PIDDV2= ( YAWE - PIDV2 ) / TT2(IT2) 

C PIDDV3= (PIDV2- PIDV3 ) / TT2(IT2) 

C D = ( KK1 ( 1 ) *KK2 ( IK2 ) *TT1 ( IT1 ) *PIDDV3 ) + PIDV4 

C C + ( KK1 ( 1 ) *KK2 ( IK2 ) *YAWE ) 

C 

c *** THE RUDDER ANGLE COULD NOT BE CHANGED MORE THAN **** 
C*** 2.5 DEGREES PER SECOND. **** 

TEM=D-TEMP 

IF (TEM. GT . 2 . 5/57 . 296) D=TEMP+2 . 5/57.296 
IF ( TEM . LT . -2.5/57.296) D=TEMP -2 . 5/57 . 296 
IF ( D.GT. (35./57.296) ) D= 35./57.296 

IF ( D.LT. (-35./S7.296) ) D = -35./S7.296 

TEMP=D 
C 

c*** THE values of hydrodynamic coefficients to represent** 

C*** SHIP STEERING DYNAMICS. ( REF. 1, 6 ) ** 

C AXIAL FORCE HYDRODYNAMIC COEFFICIENTS (SURGE) 

Xl= ( -0. 00120) *(RO*L2*S) 

X2=(0. 00045)* (RO*L2 ) 

X3= ( -0 . 000103 ) * ( RO*L2 /S ) 

X4=( -0. 00898) *(RO*L2) 
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X5= (0.00018)* (RO*L4 ) 

X6=( -0 .000948 )*(RO*L2*S*S) 

X7= ( 0.00798)*( RO*L3 ) 

X8= ( 0. 000832 )*(RO*L2*S) 

C LATERAL FORCE HYDRODYNAMIC COEFFICIENTS (SWAY) 

C Y0=( -0.000036)* (RO*L2*S*S) 

Y0=0 . 0 

Yl=( -0.011 604)*(RO*L2*S) 

Y2=( -0.08078)*(RO*L2/S) 

Y3=( -0. 000038) *(RO*L2*S) 

Y4=( -0.00499 )*(RO*L3*S) 

Y5=( 0 . 15356 )*(RO*L3/S) 

Y6= (0. 002779 ) * ( RO*L2*S*S ) 

Y7=( -0.0009)*( RO*L2* S* S ) 

Y8= (0.01 1896) *(RO*L2 ) 

C MOMENT ABOUT Z-AXIS HYDRODYNAMIC COEFFICIENTS (YAW) 
C N0= ( 0 . 000028 ) * ( RO*L3*S*S ) 

N0=0 . 0 

Nl=( -0. 002635 )*(RO*L3*S) 

N2=(0. 016361 )*(RO*L3/S) 

N3= (0. 000125 )*(RO*L3*S) 

N4=( -0.00 166) *(R0*L4*S) 

N5=( -0.05483 )*(RO*L4/S) 

. N6=( -0.001388)*( RO*L3*S*S ) 

N7= (0. 00045 )*(RO*L3*S*S) 

N8=( -0. 00489 )*(RO*L3) 

C COMMON COEFFICIENTS 

C1=(0.0084)*(RO*L3) 

C2= ( 0 . 01546 ) * ( RO*L3 ) 

C3=( -0.000086 )*(RO*L4) 

C4=(0. 000829 )*(RO*L5) 

C5= ( -0. 000227 )*(RO*L4) 

C 

C 

C REGULAR WAVES 
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FX=WA*RX*DCOS(WE*TIME+TX) 

FY=WA*RY*DCOS(WE*TIME+TY) 

MZ=WA*RZ*DCOS(WE*TIME+TZ) 



IF (TIME . EQ . 0 . 0 ) FX=0 . 0 
I F ( DABS ( FY ) . LT . 0 . 00000001 ) FY=0 . 0 
IF (DABS (MZ) . LT.O. 00000001 )MZ=0.0 



C 

c 

C EQUATIONS OF MOTION 

FI = X1*DU + X2 *DU*DU + X3*DU*DU*DU + X4*V*V + X5*R*R 
1 + X6*D*D + X7*V*R + X8* V*D + RES + FX - FXX 

F2 = YO + Y1*V + Y2*V*V*V + Y3*V*D*D + Y4*R 
1 + Y5 *R* V* V + Y6*D + Y7*D*D*D + Y8*D*V*V 

1 + FY + FYY 

F3 = NO + N1*V + N2*V*V*V + N3*V*D*D + N4*R 
1 + N5*R*V*V + N6*D + N7*D*D*D + N8*D*V*V 

1 + MZ + MZZ 



C 

C 

C*** ADDED RESISTANCE RELEVANT DUE STEERING AND TOTAL ***** 
C*** RESISTANCE AFTER THE TRANSIENT TIME. ***** 

C*** WHEN FINDING OPTIMAL PID PARAMETERS, ***** 

C*** USE (TIME. LE. 500. ) ***** 

IF (TIME. LE. 1000. ) GO TO 0001 

DELST1= ( (X4*V*V+ X5*R*R+ X6*D*D+ X7*V*R+ X8*V*D) 

1 + DELST1 ) 



TOTFR2= ( X4*V*V + X5*R*R +X1*DU +X2*DU*DU 



1 +X3 *DU*DU*DU + X6*D*D + X7*V*R + X8*V*D 

1 + FX - FXX- 63420.) + TOTFR2 



C 

c 

C*** FINDING ACCELERATION VALUES OF U, V, AND R ***** 
0001 UDOT = Fl/Cl 

VDOT = (C4*F2-C3*F3 )/(C2*C4-C5*C3 ) 



74 



RDOT 



= (C2*F3-C5*F2)/(C2*C4-C5*C3) 
C 
C 

C*** WHEN TO PRINTOUT 

IF ( I COUNT . EQ . 2 ) GO TO 50 
GO TO 300 
C 

c*** CONVERT RADIANS TO DEGREES 
50 YAWDEG= YAW*57.296 
RDEG=R*57 . 296 
RDDEG=RDOT*57 .296 
DDEG=D*57 .296 
YAWC=YAWC*57.296 
C 

C*** ONLY USE WHEN GETTING GRAPHS 
MH( ITIME+1 ) =DDEG 
T ( ITIME+1 )=TIME 
C 

ITIME=ITIME+ 1 
ICOUNT=l 
C 

C*** TEST IF WANT TO STOP 

300 IF (TIME.GT.ETIME) GO TO 400 
C 

C INTEGRATION STEP SIZE DELT 
DELT=1 . 

C 

C INTEGRATION 

u=u+udot*delt 

V=V+VDOT*DELT 

R=R+RDOT*DELT 

YAW=YAW+R*DELT 

C 

C FOLLOWINGS ARE FOR D TYPE CONTROLLER. 
S1=S1+DS1*DELT 



***** 



***** 



******* 



******* 
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S2=S2 + ( YAWE/KD ) *DELT 



C 

C 

C 

C 

C 

C 



FOLLOWINGS ARE FOR PID. 

PIDV2=PIDV2+PIDDV2*DELT 
PIDV3=PIDV3+PIDDV3*DELT 
PIDV4=PIDV4+ ( YAWE* (PIDK1/PIDT3 ) ) *DELT 



* * x * * 

* * * * * 
***** 



C*** SO MUCH SPEED DECREASING MAKE THE NONLINEAR***** 
C*** COEFFICIENT VALUES INCORRECT TO REPRESENT ***** 
c *** THE SHIP STEERING DYNAMICS. 
c *** SO RPM (PROPELLER ROTATION PER MINUTE) IS. 

C*** INCREASED TO PREVENT SHIP SPEED DECREASE. 

I F ( RPM . GT . 95 . ) GO TO 65 
I F ( U . LT . 14.7*1. 689 ) RPM=RPM+ 1./10 
65 IF(U.GT. 15.3*1. 689) RPM=RPM - 1 . /10 . 

C 

C*** TO FIND SHIP MOTIONS ON THE EARTH SURFACE. 

C*** (REFS. 1, 9) 

XDOT = U*DCOS ( YAW ) - V*DSIN(YAW) 

YDOT = U*DS IN ( YAW ) + V*DCOS(YAW) 

X = X + XDOT*DELT 
Y = Y + YDOT*DELT 



* * * * * 
* * * * * 



C***HOW MUCH ADDED FORCES OBTAINED BY INCREASING RPM***** 

RD I F=4480 . *RPM-243908 . 

RES=RD IF- 63420 . 

C*** WITH USING SHIP ENGINE SPECIFICATIONS, THE FUEL *** 
C*** CONSUMPTION CAN BE FOUND. *** 

EHP=263 . 448*RPM- 15 153 . 5312 

IF ( RPM.LE.80. ) SHP = 250*RPM - 11900. 

IF ( RPM.GT.80. ) SHP = 420*RPM - 25500. 

BARFUE = BARFUE + (0.03125*SHP +95.) * ( DELT/86400 . ) 
LBFUEL = LBFUEL + ( 0 . 000007187*SHP + 0 . 632625 ) *SHP 
C * (DELT/3600 . ) 

C 
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TIME=TIME+DELT 
I COUNT= I COUNT + 1 ' 

C*** FINDING APPROXIMATE COST FUNCTION ***** 

IF ( TIME . GT . 1000 . ) I SE=I SE + LAMDA* YAWE* *2 
I F (TIME . GT . 1000 . ) I SR=I SR + D**2 
TDI FF= ( I SE + I SR) 

C 

GO TO 200 
400 CONTINUE 
C 

c*** FINDING FUEL CONSUMPTION RATIO ***** 

DELST1 = DELST1 / (TIME - 1000.) 

TOTFR2 = TOTFR2 / (TIME - 1000.) 

PERC2 = (DELST1/TOTFR2 ) * 100. 

C 

WRITE( 6, 1111 ) TDIFF,DELST1, PERC2 
1111 FORMAT (IX , F7.3, IX, F8 . 0 , IX, F8 . 3 ) 

C 

C ***USE THIS PART TO PRINT ALL COMBINATIONS OF PID***** 

c***parameters ***** 

C WRITE (8,43) KK1 ( 1 ) , KK2 ( IK2 ) , TT1 ( IT1 ) , TT2 ( IT2 ) 

C C , TT3 ( I T3 ) , TD IFF 

C WRITE (6,43) KK1 ( 1 ) , KK2 ( IK2 ) , TT1 ( IT1 ) , TT2 ( IT2 ) 

C C , TT3 ( I T3 ) , TD I FF 

43 FORMAT (' ' , 5F8 . 2 , 2X, F20 . 5 ) 

C 

C*** DISSPLA STATEMENTS USING TO GET GRAPHS ***** 

XMIN=0 
XMAX=ETIME 
XINC=XMAX/6. 

MAGMAX= - 1 . E 1 5 
MAGMIN= 1.E15 
DO 150 1=1, IT I ME 

IF ( MH ( I ) . GT . MAGMAX ) MAGMAX=MH ( I ) 

150 IF (MH( I ) . LT.MAGMIN) MAGMIN=MH( I ) 
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CALL AREA2D ( 6 . 0 , 3 . 0 ) 

CALL XNAME( 'TIME (#SEC.&) $',100) 
CALL YNAME (' RUDDER (#DEG.&) $',100) 
CALL HEAD IN ( ' $',100,1.2,4) 



CALL HEADIN( 



SEA STATE=8 $',100,1.2,4) 



C 



CALL HEAD IN ( ' EN. ANGLE=030 SPEED=15 

D TYPE CONT. $’,100,1.2,4) 



C 



CALL HEAD IN ( 



WF = .53 ( EF =. 33 ) 
$' , 100,1.2,4) 



CALL CROSS 
CALL Y I NT AX 

CALL GRAF ( XM I N , X I NC , XMAX , MAGM I N , ' SCALE ' , MAGMAX ) 
CALL CURVE ( T , MH , I T I ME , 0 ) 

CALL RESET ( 'CROSS' ) 

CALL ENDPL(O) 

CALL DONEPL 
C 

C*** USE THIS PART TO FIND OPTIMAL PID PARAMETERS 

C IF (TDIFF. LT. SMALL ) THEN 

C SMALL=TDI FF 

C END IF 

C55 CONTINUE 

C44 CONTINUE 

C33 CONTINUE 

C22 CONTINUE 

Cll CONTINUE 

C WRITE (8,*) SMALL 



C 



C 



STOP 



END 
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APPENDIX D 



LIST OF SYMBOLS 



Commanded Heading Angle 

Yaw Angle, measured from the vertical xz plane 
to the axes of the ship; positive in the 
positive sense of rotation abiut the z-axes 



'Pe 

X 



Y 



N 



'P ' 'Pc 

Hydrodynamic force components on ship body 
( longitudinal ) 

Hydrodynamic force components on ship body 
( lateral ) 

Resultant total moments acting on a ship 
about the z-axis, yawing moment 



u 



U 1 

Au 



Velocity components of the origin of the body axes 
relative to the fluid, longitudinal component 

Commanded longitudinal velocity component 

u - u^ 

Transverse component of velocity 



r Angular velocity of yaw 

« • * 

u, v, r Acceleration components of the origin of the 

the body axes relative to the fluid (longitudinal, 
transverse, and yawing, respectively) 



<5 



Angular displacement of the rudder, measured from the 
xz -plane of the ship to the plane of the rudder 



m 



Mass of ship 
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X G 


Coordinate of the center of mass of the ship 
ship relative to body axes 


X*,Y*,N* 


• • 

Values of X, Y, and N at v=r=v=r=0 and u=u-^ 


W a 


Significant wave height 


R i 


Exciting force magnitude 


0i 


Exciting force phase angle 


w e 


Encounter frequency 


\ave 

P 


Wave circular frequency 
Encounter angle 


g 


Acceleration of gravity 


LBP 


Length between perpendiculars 


S 


(u 2 + v 2 )V2 


w 


Natural frequency of ship's steering 


OP 


Distance from the ship pivot point to the origin 
(0.3 L) 


L 


Ship length 


X 


Weighting factor 


r a' v a' 5 a 


Amplitude of r, v, and 5 


P 


Engine power 


W net 


Work done per unit mass (btu / lbm) 


m 


Fuel mass per unit time (lbm / min) 


X calm 


R = T (1 - t) , total ship resistance without 
propeller in calm water case 


t 


Thrust deduction fraction 


T 


Propeller thrust 
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2 

Partial derivative of X with respect to v 

9 

Partial derivative of Y with respect to rv 

Partial derivative of N with respect to v 

System of reference axes through origin 
of reference axes fixed in the ship 
whose direction fixed in the ship 

System of reference axes whose origin 

and direction remain fixed in the earth 

1. Signs of all directions, forces, distances, veloci- 
ties, and accelerations are positive downward along the 
z-axes, positive to starboard along the y-axes, and positive 
forward along the x-axes and similarly along the Xq yg and 

z O-axes 

2. Signs of all angles, angular velocities, angular 
accelerations, and moments are positive if clockwise when 
facing in the positive direction of appropriate axes. 






rw 

N 

vw 

x 0 ' y 0' z 0 



x, y, z 



Notes : 
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